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Decorated membrane, comprising a thin layer of elastic film with small rigid platelets fixed on top, has been 
found to be an efficient absorber of low frequency sound. In this work we consider the problem of sound 
absorption from a perspective aimed at deriving upper bounds under different scenarios, i.e., whether the sound 
is incident from one side only or from both sides, and whether there is a reflecting surface on the back side of 
the membrane. By considering the negligible thickness of the membrane, usually on the order of a fraction of 
one millimeter, we derive a relation showing that the sum of the incoming sound waves’ (complex) pressure 
amplitudes, averaged over the area of the membrane, must be equal to that of the outgoing waves. By using 
this relation, and without going to any details of the wave solutions, it is shown that the maximum absorption 
achievable from one-side incident is 50%, while the maximum absorption with a back reflecting surface can 
reach 100%. The latter was attained by the hybridized resonances. All the results are shown to be in excellent 
agreement with the experiments. This generalized perspective, when used together with the Green function 
formalism, can be useful in gaining insights and delineating the constraints on what are achievable in scatterings 
and absorption by thin film structures. 

PACS numbers: 42.25.Bs, 43.20.+g, 43.40,+s, 43.55.Ev, 46.40.-f, 62.30.+d 


I. INTRODUCTION 


A sound absorber converts the airborne acoustic energy 
into thermal motions via irreversible processes. High effi¬ 
ciency of such processes requires not only properly matched 
impedance for the absorber, but also the dissipative ability 
to absorb the incident energy. Traditional means of acous¬ 
tic absorption make use of porous and fibrous materials (T), 
gradient index materials, or perforated panels (2) with tuned 
cavity depth behind the panels. They generally result in ei¬ 
ther imperfect impedance matching to the incoming wave, 
or very bulky structures with dimensions comparable to the 
wavelength, usually on the order of meters for low frequency 
airborne sound. 

Membrane-type acoustic metamaterial, consisting of deco¬ 
rated membrane resonator (DMR) of various forms, has been 
shown to display diverse functionalities such as efficient re¬ 
flection 00, enhanced transmission |j8], reversed Doppler 
effect (9), and near-field amplification Go). The reason for 
such extraordinary behaviors can be attributed to the effec¬ 
tive negative mass density 0Q3) and/or negative refractive 
index cam introduced by the DMR’s subwavelength reso¬ 
nances. Recent works also show that the DMR can efficiently 
absorb low frequency sound even with its negligible thick¬ 
ness G3Q3, a capability clearly beyond what is achievable 
by conventional acoustic absorbers. This has been attributed 
to the high energy density of DMR’s lateral resonances. Hy¬ 
bridization of different resonances can even lead to perfect 
absorptions in the deep-subwavelength regime, thereby real¬ 
izing the exact time-reversed counterpart of an acoustic point 
source m. The latter is known to have important applica¬ 
tions for time-reversal wave technology G3. By using hy¬ 
bridized modes, an acoustic metasurface can be realized with 
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an array of such acoustic “sinks” with extraordinary sound ab¬ 
sorption characteristics m. 


In this work, we present a generalized perspective, based 
on DMR’s special geometric characteristic, for understanding 
DMR’s sound absorption behaviors. It is shown that, owing to 
membrane’s negligible thickness, there is an equality relating 
the sum of mean complex pressure amplitudes (MCPA) of the 
incoming waves, averaged over the area of DMR, to that of 
the outgoing waves. By drawing analogy to momentum con¬ 
servation law involving two equal-mass particles, it is easily 
seen that there is a component of the total incident energy, cor¬ 
responding to the center of mass motion in the two particles 
collision, which is always conserved and therefore cannot be 
dissipated. Hence only the energy in excess of the conserved 
component is available for dissipation. This observation leads 
to upper limits for DMR’s absorption performance under var¬ 
ious scenarios. In particular, for wave incident from one side 
only, only half of the incoming energy is available for dis¬ 
sipation, thus the absorption percentage cannot exceed 50% 
Da. We verify this conclusion by experimentally observing 
the sound scattering from DMRs, and the results show ex¬ 
cellent agreement with theoretical predictions. Such general 
considerations also show that the higher than 50% absorption 
can be achieved either by allowing waves to be incident from 
both sides of the membrane, such as in the coherent perfect 
absorption (CPA) scenario m, or if multiple scatterings are 
introduced into the system by introducing a back reflecting 
surface. As an example, the acoustic metasurface with hybrid 
resonances, which has been demonstrated to completely ab¬ 
sorb waves incident from one side, is analyzed in light of the 
MCPA conservation. The outcome, while matching the results 
previously reported DU. also yields additional insights. 


In what follows. Section [D] shows how the area-averaged 
response of a subwavelength membrane structure can be re¬ 
duced to a one-dimensional problem, and the related physics, 
as embedded in the Green function formalism, is consider¬ 
ably simplified as a result. In Section III we derive a general 
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FIG. 1. Schematic illustration of the scattering process from a 
DMR’s unit. With left-incident sound wave’s pressure denoted as 
p\ exp(ikoz) and that from the right as p' 2 exp(— ikoz), the scatter¬ 
ing generates two outgoing waves p\ exp(— ikoz) andp 2 exp(ikoz). 
Here ko is the sound wave’s wavevector in air, and the double 
dot-dash lines stand for periodic boundary conditions. According 
to the conservation of mean complex pressure amplitude (MCPA), 
Pi + P 2 = Pi + P 2 ■ 


conservation rule governing the scatterings from thin mem¬ 
brane structures, and show how some conclusions can be eas¬ 
ily derived from this rule without having to resort to the wave 
equation. In Section [IV] we analyze the problem of an effi¬ 
cient thin membrane absorber using the generalized perspec¬ 
tive, and verify experimentally the derived conservation rule. 
In Section[V]the problem of a membrane resonator coupled to 
a thin, sealed air cell is analyzed, leading to the condition for 
attaining perfect absorption through hybrid resonances. This 
is followed by a brief summary in Section VI that concludes 
the article. In the Appendix we derive the formulas for ob¬ 
taining the theoretically relevant parameters from four-probe 
impedance tube data. 


II. SURFACE IMPEDANCE AND DISSIPATION 


for the acoustic wave in air, where k± denotes the wavevec¬ 
tor component normal to the membrane. It follows that the 
6W component couples only to the evanescent waves, as its 
associated k± must be imaginary. In contrast, because the k|| 
components for (VF) have a distribution that peaks at kn = 0, 
it can couple to the radiation modes. Hence, if we restrict 
our considerations to only the radiation modes, then one can 
treat the problem of DMR as essentially one-dimensional in 
character. 

The surface impedance, which characterizes the DMR’s far- 
field scattering properties, can be defined by using only the 
(IV) component, given by Z = (p)/{W), with p denoting 
the total sound pressure acting on the membrane and the over¬ 
dot denoting the time derivative. We notice that, for harmoni¬ 
cally oscillating waves with angular frequency ui = 2nf, the 
time-averaged energy flux pointing into the membrane (which 
characterizes the energy dissipation rate by the membrane) 
is given by j = / Q 27r Re((p)e~ lult )Re((W)e~ ,ult )d(u;t)/(2Tr) 
m. Since the surface impedance’s real part, Re(Z), results 
from the velocity component with the same phase as pressure, 
hence this is the component that characterizes DMR’s dissipa¬ 
tive capability. It follows that j = (p) 2 Re(Z)/(2\Z\ 2 ). This 
dissipation is caused by the viscous damping in membrane’s 
relative displacements. In particular, the displacements in the 
“deaf” component 8W dominate dissipation. Therefore, there 
is a relation between R e(Z) and 8IV as shown explicitly in 
what follows. 

The surface-averaged Green function is defined by (G) = 
(W)/(p). Hence, for time-harmonic motion, we have Z = 
i/(u>(G)). For frequencies in the vicinity of one of DMR’s 
resonances, (G) is given by ET1 


(G) 


l(Wn )| 2 

Pn(^n - W 2 - 2 iwfin) ’ 


(1) 


A DMR unit comprises a uniformly stretched elastic mem¬ 
brane decorated with relatively rigid platelets. It is impor¬ 
tant for our subsequent considerations to note that the thick¬ 
ness of the membrane is negligible (~ 200 pm), while the 
cross-sectional dimension of the membrane is subwavelength 
in scale for the frequency regime of interest. Without the loss 
of generality, we can consider it as an inhomogeneous mem¬ 
brane which can vibrate and scatter the incoming plane waves 
as shown in Fig. [T] 

Central to the understanding of DMR’s acoustic behaviors 
is that only the piston-like component of its averaged displace¬ 
ment, (W), couples to the radiative wave modes. Here W 
denotes the normal displacement field of the membrane, and 
the angular bracket denotes surface averaging. However, the 
variance of the displacement, 5W = W — (IT’), is decoupled 
from the radiation modes and can be characterized as “deaf”. 
The reason of this decoupling can be seen from the Fourier 
wavevectors k|| that delineate the lateral spatial pattern of TV'. 
For the 8W, the relevant k must have magnitude that satisfy 
the inequality |k| | > 2tt/X, where A is the sound wavelength 
in air, owing to membrane’s subwavelength cross-sectional di¬ 
mension. From the displacement continuity condition and the 
wave dispersion relation, we have (/cy) 2 + ( k ±) 2 = (27r/A) 2 


where p n = f n p\W n \ 2 dV is a parameter related to the 
displacement-weighted mass density for the membrane’s nth 
eigenmode W n , p is the local mass density, Cl the volume of 
the membrane, and w n is the relevant angular eigenfrequency. 
If we denote the viscosity coefficient of the membrane as p 
ll22l . then the dissipation coefficient j3 n in Eq. (|T| is defined 
by 


/3„ = / r?(VW* • X7W n )dV/(2p n ) 

J n 

= [ p(X75W* ■ V8W n )dV/(2p n ), (2) 


In Eq. Q we have used the fact that V(W n ) = 0. It follows 
that the relevant impedance corresponding to frequencies in 
the vicinity of a resonance is given by 



Pn 

I (Wn )\ 2 



(3) 


Equations and ([3]) indicate that, although decoupled 
from radiation modes, the “deaf” component SW can still 
affect scatterings via membrane’s dissipation coefficient /?„. 
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This can be accomplished through structural designs. For in¬ 
stance, decorating membrane with semicircular platelets can 
lead to “flapping” modes having a large SW component and 
hence large dissipation m. However, a large dissipation 
sometimes may also cause impedance mismatching, thereby 
leading to reflection instead of absorption. Hence “sufficient” 
dissipation that is compatible with impedance matching is al¬ 
ways the key to an efficient absorber. 

III. A RELATION INVOLVING THE MEAN COMPLEX 

PRESSURE AMPLITUDES (MCPA) 

In this section, we will look at the problem of thin mem¬ 
brane absorption through a different perspective. As shown 
in Fig. |T] the two incoming waves counter-propagate from 
two sides, acting on the DMR with complex pressures am¬ 
plitudes p\ and p 2 . After scattering, they are converted into 
two outgoing waves with complex pressure amplitudes p° 
and p 2 . Here, the subscript “1(2)” denotes the left (right)- 
hand side region and the superscript “i(o)” stands for incom- 
ing(outgoing) waves. Notice that the sound’s energy flux in 
air, j = J 2n Re(pe~ iult )Rs(We~ iult )d(ujt)/(2 tt) = pW/2 = 
±p 2 /(2 Zq), is opposite in directions for the incident and scat¬ 
tering waves. Here Z {] denotes the characteristic impedance 
of air. Therefore, the air’s surface-averaged normal velocities 
on two sides of the membrane are given by 

(W 1 ) = 2j ± + 2 4 = ^(p\-p° 1 ), (4a) 

Pi Pi Z o 

m = 2 -4 + ^ (pii-a ), ( 4b ) 

p° 2 p' 2 Z 0 

respectively. Since the thickness of the membrane is negligi¬ 
ble, we have {W i) = (W 2 ), i.e., there is no relative motion 
between the two sides of the membrane. That immediately 
implies that the mean complex pressure amplitude (MCPA), 
p , is conserved before and after the scattering: 

P=\ [p\ + P 2 ) = \{P\+P°2)- (5) 

Conservation of the MCPA implies that the portion of in¬ 
coming energy associated with it has to be preserved from 
dissipation and transfered to radiative outgoing waves by scat¬ 
terings. This is an direct analogy to the two (equal-mass) par¬ 
ticles scattering in classical mechanics, if we regard the p’s 
as momenta of the two particles before and after the colli¬ 
sion. In that analogy Eq. ([5) expresses the fact that the cen¬ 
ter of mass momentum is always conserved before and af¬ 
ter the collision. For waves with amplitude p\ from left and 
p 2 from right, the overall incoming energy flux is given by 
[(Pi) 2 + {P 2 ) 2 }/{^Zq), while the conserved MCPA energy is 
simply j = p 2 /Zq. The component of the total energy flux 
available for dissipation, j, must be the difference between 
the two, given by 

j = (p\~p 1 2 ) 2 /(4Z 0 ). (6) 


If we denote the energy flux being absorbed as j a b, then j d b < 
j always, and the absorption coefficient, A = j a b/(j + j), has 
an upper bound of j /( j + j). 

Some easy conclusions can be stated immediately. First, 
if the incoming wave is from one side only, then p 2 = 0 and 
p = p\/2. It follows that j = j and the absorption can at most 
be 50%. Second, for A = 1, i.e., perfect absorption, one must 
have j = 0, which is only possible if p\ = —p' 2 , so thatp = 0. 
This corresponds to the CPA scenario lfl9l . In general, for a 
given ratio a = \p' 2 \/\p\ , the maximum absorption always 
occurs when p\ and p 2 are opposite in phase, so that they can 
maximally cancel each other such that 

A _ J _ {PiY + (P2) 2 + ‘APl&A 

j+j 2(^)2 +2( p ' 2 )2 

1 a 


It is interesting to see under what condition(s) would j be 
completely dissipated when the condition of a = 1 is re¬ 
laxed. That requires the total outgoing energy flux, given by 
[(Pi) 2 + (P2) 2 ]/(2 -Zo)> be equal to the conserved MCPA en- 
ergy j = p 2 /Zq. This can occur when p° = p 2 = p. Since the 
two outgoing waves are opposite in direction, their pressure 
exactly cancel each other on the membrane. Hence the net 
pressure on the membrane is given by (ptot) = p\ — P 2 , while 
displacement continuity means that the membrane’s surface- 
averaged normal velocity ( W) must be equal to that of air, 
{ W \), which is defined by Eq. ( |4a| >. Thus ( W ) = ( Wi) = 
{Pi-pD/Zo = {p\-p l 2 )/{2Z 0 ) since p\ =p = {p\+pl 2 )/2. 
In other words, the requirement boils down to 

z = (Ptot)/(WO = 2 z 0 . (8) 

Therefore, as long as a DMR has a purely real surface 
impedance twice the impedance of air, it has the ability to 
absorb all the available energy j in the incoming waves. Ac¬ 
cording to Eq. 0- such impedance is realizable at the reso¬ 
nances of a DMR, i.e., when u> = u> n , with a suitably valued 

Pn- 

IV. DECORATED MEMBRANE WITH SEMICIRCULAR 
PLATELETS 

In this and the following Sections, we shall present two sce¬ 
narios involving large dissipation by the DMRs in which the 
absorption measurements are performed by using the setup 
shown in Fig.[2ja). In the present Section we consider absorp¬ 
tion by DMRs with asymmetric platelets. 

As the rubber membrane is usually weakly dissipative in 
character, we decorate it with semicircular platelets to utilize 
the high concentration of curvature energy at the edges of the 
platelets to enhance the dissipation. The asymmetry of the 
platelets means that there can be “flapping” modes [Figs.[2jc) 
and (d)] which represent a combination of normal displace¬ 
ment of the platelets in conjunction with a rotational motion. 
At the perimeters of the platelets, there can be very high con¬ 
centration of so-called curvature energy that is proportional to 
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FIG. 2. (a) Schematic illustration of the experiment apparatus. The 
front tube with a loudspeaker at one end has two sensors, while the 
back tube having the other two sensors is plugged by acoustic foam 
at another end. These four sensors are marked from 1 to 4 in suc¬ 
cession. (b-e) Front views of the DMR's structure and its eigen- 
modes’ profiles. Decorated by two opposite semicircular platelets 
[gray semicircles in (b)], the stretched membrane [white rectangle in 
(b)] presents three eigenmodes below 1 kHz whose normal displace¬ 
ment amplitude | | is plotted in gray scale as shown in (c-e). 



FIG. 3. Evidence of MCPA conservation and total absorption of 
available energy under one-side incident wave, (a) Small discrepan¬ 
cies between R + T and 1 (< 2.97%) confirm the conservation of 
MCPA. Equal amplitude, close to 0.5 for both T in (b) and R in (c) 
(only their magnitudes have been shown here for clarity) at 383 Hz, 
indicates all the available energy has been absorbed, as shown by the 
absorption peak of 49.327% at the same frequency in (d). The red 
circles are the experimental results, the black curves are calculated 
from the eigenmodes shown in Figs.|2|c-e). The inset shows the ex¬ 
perimental sample with the red square indicating the area covered by 
the impedance tube. 


the square of the second derivative of the normal displacement 
along the membrane surface directions. Since the energy dis¬ 
sipation is given by the integral of the product of energy den¬ 
sity with the dissipation coefficient, a very high energy density 
can significantly enhance the absorption of sound. Moreover, 
with platelets of different weights, the size of this edge effect 
can be tuned to give different R e(Z). 

Fig. |2|b) shows a unit of our DMR. Four edges of the 
stretched rectangle membrane are fixed on a rigid frame. On 
top of the membrane, two decorated semicircular platelets are 
fixed to face each other. The membrane has a width of 15 mm, 
length of 30 mm, and a thickness of 0.2 mm. While the mass 
of platelets can vary in different samples, the radius for each 
platelet is kept at 6 mm. 

We test the DMR’s scattering and absorption properties by 
mounting a unit of the DMR on a rigid plate, and the whole 
sample is sandwiched by two impedance tubes having square 
cross-sections (see the inset in Fig. [3J. The front tube has 
two sensors [Fig. [2ja)], plus a loudspeaker at the front end 
to generate the plane waves. The back tube has another two 
sensors, and the tube’s back end is filled by acoustic foam to 
eliminate reflection. By normalizing the pressure amplitude 
of all the relevant sound waves by the incident sound pressure 
amplitude, the reflection and transmission coefficients, R and 
T, can be obtained from the pressure data recorded by the four 
sensors (see Eq. |A5| in Appendix). 

In the present case the incident wave is from one side only, 
hence in accordance to their definitions, R = p\/p\ and 
T = j >2 /p\ ■ The MCPA conservation in this case is given 
by the simple relation 1 = R + T (notice that this is not the 


energy conservation law, which would involve the square of 
R and T). We can check it directly by comparing the mea¬ 
sured values of R + T with 1. For a sample with identical 
platelets’ mass of 151 mg each. Figs. [3ja-c) show the mea¬ 
sured R + T, |T|, and \R\ respectively. The small discrepancy 
between R + T and 1 (< 2.97%) confirms the validity of the 
MCPA’s conservation over all frequencies of interest. 

We notice that, at the resonant frequency of 383 Hz, both R 
and T have an approximately equal value (~ 0.5). It indicates 
an almost total absorption for the available energy j (which 
equals to the MCPA energy j here), hence the absorption co¬ 
efficient A reaches its maximum value ,4 max =0.5 [i.e., a = 0 
in Eq. 01- This is seen from the measured absorption coeffi¬ 
cient A = 1 — \R\ 2 — |T| 2 with a peak magnitude of ~ 0.493 
at the same frequency as shown in Fig. [3jd). 

The displacement profile of this (maximum absorption) res¬ 
onant mode is shown in Fig.[2jc) [together with Figs.[2[d) and 
(e) for the other two absorption peaks in Fig. [3}d)], which 
is numerically simulated (based on the material parameters 
in Ref. f2T|) by using COMSOL Multiphysics—a commer¬ 
cial finite-element solver software. Based on this eigenmode, 
the evaluated parameter values from the experimental data 
are \(Wi)\ 2 / Pl = 0.039m 2 /kg, wi = 2 tt x 283.43 Hz, 
pi = 15.74 Hz and // = 1.76 Pa-s. The surface impedance of 
this DMR can therefore be calculated by using Eq. |3]i to yield 
Z = 1.986%), which is very coincident with the condition as 
expressed by Eq. ([8]). Resonances’ departures from this condi¬ 
tion are caused by mismatching in impedance, thereby leading 
to weaker absorptions. We demonstrate this effect by 7 sam¬ 
ples having platelets’ mass varying from 34 to 825 mg, and 
comparing the peak absorption coefficient ,4 pcak with the rel- 
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FIG. 4. Maximum absorption for impedance-matching and its re¬ 
lation with incident waves’ symmetry, (a) To absorb one-side inci¬ 
dent wave, the DMR samples with differently weighted platelets are 
shown to display various peak absorption coefficients Apeak, with the 
maximum achieved when the relevant impedance satisfies the condi¬ 
tion as expressed by 0 — Z = 2Zq, with Z Q being the characteristic 
impedance of air. The black lines are from numerical simulations 
while the red circles denote experimental data, (b) Different max¬ 
imum absorption coefficient A milx for two-side incident waves for 
various values of a 2 . The numerically simulated result (red circles) 
matches Eq. 0 (black curve) very well. 


evant surface impedance in Fig. |4ja). Clearly, the maximum 
peak absorption is attained when the condition expressed by 
([8J is met. 

Through numerical simulations, we further locate the max¬ 
imum absorption coefficient A max for two-side incoming 
waves with various values of a 2 = [pil' 2 /A 2 > by using 
the DMR with platelets’ mass of 151 mg each. As shown 
in Fig. |4jb), almost perfect agreement between the numeri¬ 
cally simulated results and that predicted by Eq. 0 is seen. 
This confirms again the validity of the “available energy for 
dissipation” concept, and its relation to MCPA. 

Combinations of different units in one panel can give rise 
to a broader absorption spectrum. We demonstrate this by us¬ 
ing a panel consisting of eight different units (see the inset of 
Fig. 0- Five of them have decorated platelets with mass of 
225 mg each, and the other three have platelets with masses 
of 70, 140, and 445 mg. The experimentally measured ab¬ 
sorption coefficient A for this panel is shown in Fig. [5] A 
relatively broad absorption spectrum [compared to that shown 
in Fig. [3jd)] is seen due to the merging of multiple absorption 
peaks. We would like to note here that none of the absorption 
spectra in Fig. [3jd) and Fig.[5]exceed 50%. In this context we 
would like to note that in Ref. era, maximum absorption over 
70% has been reported for a one-side incident wave configura¬ 
tion. The reason for this discrepancy is simply the absorption 
area. The DMR panel in Ref. Q4] has a larger area (about 
2.26 times larger) than the incoming wave front, and part of 
the sample extended outside of the impedance tube. Here, 
the tested samples are always contained within the impedance 
tube, and have exactly the same cross-section [see the insets 
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FIG. 5. Experimentally measured absorption spectrum for a DMR 
panel with 8 units. Five in these eight units have the same deco¬ 
rated semicircular platelets weighting 225 mg, while the other three 
platelets have masses of 70, 140, and 445 mg. The inset shows the 
experimental sample, in which the red square indicates the area cov¬ 
ered by the impedance tube. 
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FIG. 6. Schematic illustration of the multiply-scattered waves 
that can lead to perfect absorption on a DMR, with a reflecting 
wall placed a small distance s behind it. (a) When the incom¬ 
ing wave energy has been totally absorbed, there should be no re¬ 
flections, i.e., all the scatterings on the incident side cancel each 
other pi exp (—ikoz) = X^m=o Pi 7 " exp(— iko z ) = 0. Incom¬ 
ing wave pi exp(ifcoz) generates only two net counter-propagating 
wave components p 2 exp(ifco-z) = Yim-oPi"" ex P(^oz) and 
p 2 exp(— ikoz) = Ylm=l Pi"' exp(— ikoz), in the air layer between 
the DMR (at z = 0) and the reflecting wall. Here, the double 
dot-dash lines stand for the location at which the periodic boundary 
condition is applied, (b) Due to the MCPA conservation, the three 
complex pressure amplitudes must form an isosceles triangle whose 
two equal sides are slanted toward each other with an angle given by 
25 = 2kos. 


in Fig. [3jd) and Fig. [5]. 


V. PERFECT ABSORPTION BY ACOUSTIC 
METASURFACE 


Although the MCPA energy in the incoming waves has to 
be conserved in a single scattering event, its dissipation is still 
possible through multiple scatterings. The outgoing MCPA 
energy after each scattering can serve as the incoming energy 
again for the subsequent scatterings. An example is the acous¬ 
tic metasurface reported in Ref. ed. in which a reflecting 
wall is placed behind a DMR, separated by a distance s that 
is deeply subwavelength (2 orders smaller than the relevant 
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wavelength in air), thereby creating multiple reflections be¬ 
tween the membrane and the reflecting wall [see in Fig.JhJ'a)]. 
As a result, the absorption coefficient for one-side incident 
wave can reach nearly 100%. In this section, we analyze this 
multiple-scattering process and the perfect absorption condi¬ 
tion with the aid of MCPA conservation. 

In the air layer between the membrane and the reflecting 
surface, the outgoing wave after the mth scattering from the 
membrane, p 2 m , is reflected by the wall, and becomes the in¬ 
coming wave, p 2 m+1 > for the (to + l)th scattering. These mul¬ 
tiple scattered waves can be superposed to form two counter- 
propagating waves: p\ = £~ =1 'PT and p° = Y,Z=oP° 2 m - 
Because the reflecting wall presents a velocity node, these two 
wave components in front of the reflecting wall must be equal 
in magnitude but different in phase by 25 = 2k^s at the posi¬ 
tion of membrane, i.e., 

\P 2 1 = \P° 2 l P 2 /P 2 = exp(2 iS). 

Here ko denotes the sound’s wavevector in air, and s is the air 
layer’s thickness. 

As the conservation of MCPA is valid for membrane’s each 
scattering, it also holds for their superpositions—the total 
wave field on the membrane. Therefore, (p\ + p‘ 2 )/2 = 
{p° + P° 2 )/2, where p\ = Em=oPi m * s a ^ so a superposition 
for outgoing waves on the left-hand side after each scatter¬ 
ing. Assuming the incoming energy is completely absorbed 
by the DMR, there will be no reflection on the incident side, 
thus p° = 0. Thus the conservation of MCPA, p\ = p 2 - pl 2 , 
may be represented by a phasor-diagram, shown in Fig. [6j b), 
which is an isosceles triangle with p l 2 and p 2 being the two 
equal sides and p\ being the base. It is clear from the diagram 
that p 2 +p 2 , which is the right-hand side pressure on the DMR, 
differs from the left-hand side pressure, p\ , by a phase of 7r/2. 
That is, p 2 + p 2 is perpendicular to p\ in the phasor-diagram, 
so that p 2 + p 2 = ip\ cot 5. Since the surface-averaged ve¬ 
locity is given by (W) = (W\) = p\/Z 0 , the total pressure 
applied to the DMR, (p tot ) = p\ — (p 2 + p 2 ), is given by 
(Ptot) = p\ (1 — i cot 5). It follows that the impedance condi¬ 
tion for achieving total absorption is given by 

Z = (p M )/(W) = Z 0 (l-icatS). (9) 

For 5 = 7t/2, we have Z = Zq, i.e., impedance matching. 
Such condition occurs at the sealed cell’s “drum” resonances 
when s = (N/2 + 1/4) A, with N being an integer. However, 
for an air layer thinner than a quarter-wavelength, <5 —* 0, and 
the required imaginary part of Z approaches —ioo. Such a 
large imaginary part of the DMR impedance is necessary in 
order to cancel the impedance of the reflecting back surface 
that is close by. Since the impedance of the whole structure 
is the addition of the two serial impedances of the DMR and 
the sealed cell behind it, the net impedance of the structure is 
still Z = Z 0 , i.e., impedance matched to that of air fi~8l . In 
what follows, we will show that by utilizing the hybrid mode 
of the DMR, such large imaginary impedance of the DMR can 
indeed be realized, while maintaining the real part to be Zq. 
The reason behind the large imaginary impedance is due to 
the existence of DMR’s anti-resonance condition, realized at 




FIG. 7. The behavior of a perfect absorber for one-side incident 
wave, denoted the hybrid resonance metasurface, (a) Schematic il¬ 
lustration of the unit cell’s component and geometry, (b) Schematic 
cross-sectional illustration of the two lowest frequency eigenmodes 
of the DMR. plotted in reference to the same phase of incident wave, 
with W being the normal displacement of membrane, normalized to 
its surface average component (W). (c) Displacement profile of the 
DMR’s hybrid mode at 152 Hz, which is clearly a superposition of 
the two in (b), but with an almost an order of magnitude larger ampli¬ 
tude. However, such large amplitude oscillations are only coupled to 
the non-radiative evanescent waves, as shown in Sectionjll] (d,e) As a 
function of frequency, the DMR’s surface impedance Z , normalized 
to the impedance of air, Zo, is seen to satisfy the total absorption re¬ 
quirement given by Eq. at the hybrid resonance frequency, which 
is indicated by the blue dot-dashed line. This is confirmed experi¬ 
mentally measuring the absorption coefficient, shown in (f). A sharp 
absorption peak, reaching 0.994, is seen. The solid and dashed ar¬ 
rows indicate the first eigenmode and the anti-resonance frequencies 
of the DMR. Experimental results are shown in red circles, and black 
curves denote the simulation results. 


frequency u) that is in-between the two low-lying resonances, 
at which the dynamic mass density of the DMR can display a 
resonance-like dispersion 0. 

In order to see the emergence of the hybrid mode, here 
for simplicity we consider only two relevant eigenmodes IT-| 
and IV 2 [as shown in Fig.[7Jb)], in contrast to Eq. (|T|, which 
is valid only for frequencies close to one of the DMR reso¬ 
nances. The surface-averaged Green function of the DMR is 
given by J2T1 


2 


<G) = E 

n—1 


l(^n)| 2 

Pn{ul - W 2 ) 


2 

+2/3 

n—l 


\(W n )\ 2 u 

Pn(u* - W 2 ) 2 


( 10 ) 


where /3, the averaged coefficient of B\ and B 2 , is taken to be a 
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small quantity so that the dimensionless P/u <C 1. We notice 
that in the vicinity of the DMR’s anti-resonance frequency A, 
at which Re((G)) = 0 (i.e., Z —> oo), this surface-averaged 
Green function behaves as ( G) ~ 2S(*/3 — Aw), where 


|(W„)| 2 w 


and Aw = w — w. 


From the inverse relationship between the Green function and 
impedance as describe previously, the surface impedance in 
the vicinity of DMR’s anti-resonance is given by 


Z = 


1 1 

2wS f3 2 + Aw 2 


(/3 — iAw). 


( 11 ) 


Equating Eqs. |9]l and ( fTT) yields two equations for the condi¬ 
tion of total absorption: 


1 P 

2wS /3 2 + Aw 2 


— Z 0 , 


1 Aw 
2wS /3 2 + Aw 2 


Z 0 cot S. 


(12a) 

(12b) 


This total absorption condition is robust since experimentally 
we have two parameters, fj and Aw, which can be adjusted 
to satisfy Eq. Notice that this total absorption mode of 
the membrane is neither of its two original DMR resonances, 
but a hybridization of them [as shown in Fig. [7jc)] caused by 
the back reflecting surface that is located in the near-field re¬ 
gion of the DMR. Since the dissipation coefficient is required 
to be very small, hence to completely dissipate the incoming 
acoustic energy, the amplitude of this hybrid mode is always 
significantly larger than the original membrane’s resonances 
for the same incident wave amplitude. This is clear from the 
comparison between Figs.[7Jb) and (c). However, it is to be 
especially noted that in spite of the large maximum amplitude 
of the hybrid mode, the surface-averaged normal displacement 
is small and matches that of air. 

Experiments reported in Ref. Ifl8l have corroborated the 
existence and condition for this hybrid resonance total ab¬ 
sorption. By placing an aluminum-reflecting wall behind 
a DMR, which is a 90 mm-wide, tensioned circular elastic 
membrane decorated by a 800 mg platelet with a radius of 10 
mm [Fig. [7ja)], an absorption peak with A > 0.99 was ob¬ 
served at 152 Hz [Fig. [7}f)]. By using a value of s = 21.7 
mm, which is inferred from the experimental data of s = 17 
mm of SF6 gas, as the effect of two thin gas layers is iden¬ 
tical if their thicknesses are scaled linearly with their adia¬ 
batic index (SF6 has an adiabatic index of 1.098, compared 
to 1.4 for air) ED- Based on the two relevant eigenmodes 
shown in Fig. |7jb), the evaluated S = 1.268 x 10 _9 m 2 s 3 /kg, 
w = 2-7T x 162.3 Hz, and /3 = 5.14 Hz. Equation © gives 
the surface impedance of the DMR as a function of frequency, 
shown in Figs. [TJd) and [7J e). At the total absorption fre¬ 
quency, we have Z = 1.001 — 12.961*, which is noted to 
match the value required, Z 0 (l — icot(fcos)) = 1 — 16.342i, 
reasonably well. 


VI. CONCLUDING REMARKS 

We have proposed a generalized perspective for the scatter¬ 
ing and absorption behaviors of thin membrane structures in 
the subwavelength regime. As a result of DMR’s negligible 
thickness, the conservation of MCPA for the incoming and 
outgoing waves can be derived. By analogy to the momen¬ 
tum conservation law involving two equal-mass particles, it is 
easy to see that within one scattering, waves incident from one 
side can only be absorbed by a maximum of 50%. Higher ab¬ 
sorption requires waves incident from both sides of the DMR. 
Based on the Green function formalism, we determined the 
value of maximal absorption under different scenarios. In par¬ 
ticular, maximal absorption is realized when the DMR has a 
surface impedance twice that of air. These conclusions are ex¬ 
amined experimentally. In all cases good agreement with the 
theory is obtained. 

For one-side incident, absorption higher than 50% can be 
achieved by introducing multiple scatterings, e.g., by plac¬ 
ing a reflecting surface behind the DMR. Analysis by using 
the MCPA conservation law shows that near-perfect absorp¬ 
tion can be achieved through hybrid resonances, in excel¬ 
lent agreement with the experimental results reported recently 
ED; while at the same time it also yields the general condi¬ 
tion [Eq |9]l] for achieving such characteristic. 

Due to the similarity between the acoustic and electromag¬ 
netic waves, there can clearly be analogous considerations for 
the latter, e.g., in laser plasma absorption |[23l425l . Hence sim¬ 
ilar results are expected for electromagnetic wave scatterings 
from thin film structures. In particular, similarity may extend 
to perfect absorption by electromagnetic hybrid resonances. 
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Appendix A: Transmission/Reflection Retrieval Method 

As scalar waves, airborne sound can propagate in a sub¬ 
wavelength waveguide without a cut-off frequency. In our 
experiments, the geometrical size of the apparatus (viz., the 
width of the square waveguide) [cf. Fig.[2](a)] is smaller than 
the measured wavelength, so that only plane waves can prop¬ 
agate in both the front and back tubes (26). The total pressure 
fields in the two (front and back) impedance tubes may be 
expressed as the sum of forward and backward waves propa¬ 
gating along the z direction: 

Pl = p\e ikoZ + ple-^, (Ala) 

p 2 = p' 2 e~ ikoZ + p° 2 e ikoZ . (Alb) 
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Here the subscripts “1” and “2” refer to the front and back 
tubes and the superscripts “i” and “o” represent the incoming 
and outgoing waves, respectively. To retrieve the transmission 
and reflection coefficients, the total pressure fields should be 
exactly expressed, by using the four experimentally measured 
parameters: p \, p °, p \ 2 and p 2 . Four sensors are used to deter¬ 
mine these parameters. Two sensors labeled “1” and “2” are 
placed in the front tube at Z\ = —339.5 mm and z 2 = —239.5 
mm, and the other two sensors labeled “3” and “4” are placed 
in the back tube at z 3 = 193.0 mm and Z 4 = 393.0 mm [cf. 
Fig. [2] (a)]. According to Eq. ( |A1| >, the pressure values at the 
positions of the four sensors are 


P(z 1 ) =P\e 
p{z 2 ) = p\e 
p{z 3 ) = P 2 e 
p{z 4 ) = P 2 e 


ik 0 z 1 _|_ p^ e ~ ik o z i 
ik n zn 1 o —iknZn 

u + p 3 e u 
-ik 0 z 3 -p p° e ik 0 z 3 

-ik 0 z4 pO e ik 0 z 4 


By solving Eq. (|A2[>, we obtain 


P\ 

Pi 


P\ 

P °2 


p(zi)e ik ° Zl - p(z 2 )e ik ° z2 

g2ikoZi _ g2ikoZ2 

p(zi)e ik ° Z2 — p(z 2 )e ik ° Zl c ^. n ( Zl+Z! ,) 

g2ikoZi _ g2ikoZ2 

p(z 3 )e ik ° Z4 -p(z 4 )e lk ° Z3 c , fcn(z , +Z4) 

g2ikoZ3 _ g2ikoZ4 

p{z 3 )e lk ° Z3 — p(z4)e lk ° Zi 

g2ikaZ3 _ g2ikoZ4 


(A2a) 

(A2b) 

(A2c) 

(A2d) 


(A3a) 

(A3b) 

(A3c) 

(A3d) 


Here p(zj ) is the pressure measured by each sensor labeled as 
“j = 1 ~ 4” in the subscripts. 

The scattering matrix S(ko) describing the relationship be¬ 
tween the incoming and outgoing waves can be expressed as 


S{ko) 


T R 
R T 


(A4) 


It should be noted that, due to the symmetry of the sample 
in our system, the reflection and transmission coefficients R 
and T are identical if the sample is tuned around 180 degrees. 
As a result, the reflection and transmission coefficients can be 
retrieved as 


R = 

T = 


P1P1 - P2P2 

P1P1 - P2P2 ’ 

(A5a) 

P1P2 - P1P2 

P1P1 - P2P2 ' 

(A5b) 


Here the four wave amplitude p\, p '\, p ' 2 and p 2 are determined 
from Eq. ( |A3| ). 
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